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ABSTRACT 

v 

'  This  paper  considers  the  rescheduling  of  operations  with  release  dates  and  multiple  resources 
when  disruptions  prevent  the  use  of  a  preplanned  schedule.  The  overall  strategy  is  to  follow  the 
preschedule  until  a  disruption  occurs.  After  a  disruption,  part  of  the  schedule  is  reconstructed 
to  matchup  with  the  preschedule  at  some  future  time.  Conditions  are  given  for  the  optimality 
of  this  approach.  A  practical  implementation  is  compared  with  the  alternatives  of  preplanned 
static  scheduling  and  myopic  dynamic  scheduling.  A  set  of  practical  test  problems  demonstrates 
the  advantages  of  the  matchup  approach.  We  also  explore  the  solution  of  the  matchup  scheduling 
problem  and  show  the  advantages  of  an  integer  programming  approach  for  allocating  resources  to 
jobs. 

Keywords:  scheduling,  integer  programming,  unreliable  machines. 

1.  Introduction 

Much  of  the  research  in  scheduling  considers  environments  with  one  or  more  of  the  follow¬ 
ing  assumptions:  a  single  required  resource  (a  machine),  identical  resources,  all  jobs  available  at 
the  same  fixed  time,  or  known  and  fixed  future  conditions  (see,  e.g.,  Graves  [1081]  and  Rinnoov 
Kan[1076]).  These  conditions  are  seldom  true  in  actual  production  facilities.  In  this  paper  we 
consider  a  method  for  adapting  a  preplanned  schedule  to  a  changing  scheduling  environment.  The 
problem  includes  multiple  resources  with  some  degree  of  processing  compatibility,  varying  times 
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at  which  jobs  are  available  and  costs  associated  with  not  completing  jobs  by  a  due  date.  This 
situation  models  two  actual  manufacturing  facilities  in  the  automobile  industry  that  are  used  in 
testing  the  approach. 

Few  previous  studies  have  considered  the  possibility  for  disruptions  such  as  machine  break¬ 
downs.  When  disruptions  are  considered  (see  Glazebrook  [1084] ,  e.g.),  the  models  have  limiting 
conditions.  Models  without  disruptions  include  using  a  common  deadline  (Root  [1065]),  identical 
processors  (Dogramaci  and  Surkis  [1070]),  and  unit  processing  times  (Blazewicz  [1070]).  The  ob¬ 
jectives  studied  include  minimizing  maximum  completion  time  (Bratley,  et  al  [1075],  Nichols,  et 
al.  [1078])  or  maximum  tardiness  (Nunnikhoven  and  Emmons  [1077]). 

We  model  a  general  multiple  resources  system  with  disruptions  and  assume  that  the  presched¬ 
ule  can  be  followed  if  no  disruptions  occur.  The  scheduling  strategy  in  this  paper  follows  this 
previously  planned  schedule  until  a  disruption  occurs  and  then  reschedules  part  of  the  preschedule 
to  accomodate  the  disruption.  We  reschedule  to  match  up  with  the  preschedule  at  some  time  in 
the  future,  that  is,  to  reschedule  so  that  the  completed  jobs  and  inventory  positions  are  identical  to 
what  would  have  occurred  in  the  preschedule.  Were  the  problem  formulated  as  a  dynamic  program, 
the  state  reached  by  the  revised  schedule  is  the  same  as  that  reached  by  the  original  schedule. 

The  matchup  method  is  similar  to  the  approaches  in  Chang,  et  al.  [1084],  Donath  and  Graves 
[1085]  and  Filip,  et  al.  [1083],  in  that  existing  schedules  are  revised  at  disruptions.  Our  matchup 
method  differs  from  these  approaches  in  that  it  seeks  to  matchup  with  the  preschedule. 

We  use  a  set  of  real  problems  to  demonstrate  the  effectiveness  of  this  procedure  in  comparison 
with  pushing  back  a  preplanned  schedule  and  fully  dynamic  scheduling.  We  also  present  heuristic 
procedures  for  solving  the  matchup  problem  and  compare  them  on  the  problem  set. 

The  general  matchup  method  is  described  in  Section  2,  where  we  show  the  optimality  of 
the  matchup  approach  for  sufficiently  spaced  disruptions.  Techniques  for  sequencing  jobs  on  one 
resource  appear  in  Section  3.  Methods  for  allocating  jobs  to  the  resources  appear  in  Section  4. 
Section  5  describes  the  experiments  conducted  on  our  practical  test  set.  Section  6  summarizes  our 
results  and  conclusions. 

2.  Matchup  Scheduling  Strategy 

In  the  extremes,  scheduling  systems  may  be  run  continuously,  as  in  the  dynamic  algorithm>  nt 
Baker  and  Kanet  [1083]  and  Morton  and  Rachamadugu  [1083],  or  very  infrequently  os  in  the 
algorithms  of  Lageweg,  et  al.  [1077]  and  Bratley,  et  al.  The  matchup  method  fits  between  t 
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common  approaches.  It  responds  to  disruptions  as  in  a  dynamic  algorithm,  yet  considers  future 
information  as  in  a  static  algorithm.  Under  certain  conditions  the  matchup  approach  leads  to  an 
optimal  schedule  following  a  disruption. 


I 


I 


The  theoretical  basis  for  the  matchup  strategy  is  an  extension  of  economic  turnpike  results 
(see  McKenzie  [1976]  ).  In  Bean  and  Birge  [1985],  it  was  shown  that  a  scheduling  model  could  fit 
the  framework  of  McKenzie’s  model  for  the  asymptotic  stability  of  optimal  solutions.  We  extend 
these  results  below  by  providing  a  general  scheduling  model  that  yields  asymptotic  stability  results 
under  fairly  general  conditions. 

Most  scheduling  models  assume  a  finite  number  of  jobs  in  the  system.  However,  in  most  real 
problems,  as  jobs  are  processed  other  jobs  are  introduced.  Failure  to  recognize  the  ongoing  nature 
of  the  problem  constitutes  a  significant  simplification.  To  model  the  indefinite  time  horizon  of 
realistic  scheduling  problems  we  consider  a  discrete  time  infinite  horizon  optimization  problem.  At 
time  t  =  0, 1, . . . ,  we  are  in  state  xt  €  SR".  In  the  scheduling  framework  a  state  includes  information 
on  the  status  of  each  job,  e.g.  ,  what  fraction  has  been  completed  and  what  resources  are  assigned. 
The  sequence  of  states  over  time  is  i  €  £oc.  The  decision  at  time  t  is  to  choose  a  state  xe+l  to 
enter  at  time  t  +  1  with  cost  /t(xt,xt+j).  The  objective  is  to  find 

PC 

mf  /«(*«, Xf+i ),  (P) 

tet 

where  each  ft  is  a  proper  convex  function  of  xt  and  x(+i- 

We  assume  that  each  component  i  of  it  corresponds  to  production  of  part  i.  There  are  F 
types  of  resources  potentially  required  to  produce  each  part.  Each  resource  set  r  has  A/(r)  groups. 

G J, . . .  of  identical  resources.  Part  t  requires  a  resource  in  group  Grg (|  r)  for  r  =  1 . R.  for 

processing. 

Each  part  i  also  has  an  associated  processing  requirement  p(i),  an  earliness  cost  weight,  u(i)  > 

0.  and  a  lateness  cost  weight,  w(i)  >0.  There  is  also  a  sequence  of  release  dates,  {r(t,  1 ).  r(i,  2) _ }. 

early  dates,  {«(*,  l),e(i,2),. . .},  and  due  dates,  {d(i,  l),<f(t,2),.. .},  such  that  r(i,k)  <  e(i,k )  < 

d{i,k)  for  each  k  =  1,2 .  The  objective  function  /  has  two  components,  /,  and  6r.  The  f, 

component  is  separable  in  i  =  1,2, . . . ,n.  It  includes  a  lateness  penalty  if  the  kth  shipment  of  a 
part  is  not  completed  by  the  kth  due  date.  It  includes  an  earliness  penalty  (holding  cost)  if  the 
kth  shipment  of  a  part  is  completed  before  the  kth  early  date.  The  6r  component  forces  feasibility 
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of  processing  with  respect  to  the  finite  resources.  The  objective  function  is  then 


n  A 

ft(xt,Xt+i)  =  £/7(x, (>),*<+!(*))  + 

i=l  r=  1 


where  d{i,l  —  !)<<<  d(i,l),  and 


/i(xi(*)»*«+i(0) 


'  -WiXt(l),  if  —  l[t=<i(i,/)]P(l)  <  x<+i(t)  -  x,(i)  <  l[t>r(,,0]  -  = 

xt(i)  <  0 

l[««(i./)]u»(x‘(*)  “  P(*))i  -f[t=d(i,0]P(0  -  xt+i(*)  ~  x((i)  <  l[«>r(.,/)]  -  l[/=^i 

<  xt(i)  >  p(i) 

if  —  P(*)  <  xt+l(0  -  X«(l)  <  l(t>  r( .,/)]  -  1[(=</| 

o  <  Xt(i)  <  p(i) 

.  oo,  otherwise 


Note  that  f\  is  a  piecewise,  linear  convex  function  (depending  on  f)  of  x,(i)  for  each  feasible 
x(+i(i)  and  that  the  set  of  feasible  x<+i(i)  is  convex.  Hence  f\  is  convex.  Note  also  that  /t*  forces 
a  shipment  on  each  due  date  (possibly  resulting  in  negative  inventory)  so  that  feasible  xt  remain 
bounded.  It  represents  incremental  penalties  for  tardiness  and  earliness.  The  resource  constraints 
are  represented  by 


6r(lt,X|+i)  =  /  if  l[(K».r)=>](x«+l  x«)  +  p(i)f[<=i(i,/)]  —  l^jl>  j  —  1 . *^<r) 

\  1,  otherwise 

These  constraints  limit  the  number  of  resources  used  of  each  type  to  the  cardinality  of  those  sets. 
Since  each  of  these  indicator  functions  is  convex,  ft  is  convex. 

Problem  P  is  then  a  convex  optimization  problem.  The  objective  may  not,  however,  be  finite 
for  any  x  €  i°°.  We  avoid  this  difficulty  by  defining  a  policy  x*  =  {xo,xJ, . . .}  as  weakly  optimal,  as 

in  McKenzie,  if  x*  is  not  overtaken  by  any  other  policy,  i.e.,  if  there  does  not  exist  x1  =  (x0 , x J _ } 

such  that 

r 

limsup  y^(/«(x(,xj+1)  -  /i(xj,x;+1))  <  -e,  (I  j 

r-°c  (.0 

for  some  «  >  0.  The  following  well  known  technique  gives  us  a  finite  valued  problem.  By  the 
construction  of  ft  we  know  that  the  objective  terms  are  bounded  for  each  feasible  (x(. x,+  1  ).\Ve 
can  formulate  an  equivalent  problem  by  subtracting  the  period  t  value  of  the  optimal  solution  from 
all  ft  in  each  period.  This  altered  problem  has  the  same  optimal  solutions  as  the  original  problem 
and  is  finitely  valued  for  any  weakly  optimal  policy.  Therefore,  we  can  assume  without  loss  of 
generality  that  P  has  a  finite  optimum  obtained  by  rgf*. 


To  obtain  results  on  the  stability  of  solutions,  we  first  present  conditions  for  a  supporting  price 
system.  For  this  development,  we  define  the  set  of  feasible  states  at  time  t,  .Yt,  as  all  states  which 
can  be  continued  over  the  infinite  horizon  at  finite  cost. 

The  attainable  states  at  t  are 

Yt  =  {xtlft-ilxt-iyXt)  <  oo  for  somex<_i}. 

A  resource  group,  Gj ,  is  idle  at  t  if  no  processing  occurs  for  parts  requiring  G r.,  i.e.,  if  x(+i(e)  = 
x<(t)  -  l[t=<t(i./)]P(t)  for  all  g{i,r)  =  j. 

The  next  theorem  expands  on  McKenzie’s  results  for  price  supports  by  using  specific  properties 
of  the  objective  functions  /<. 

Theorem  1:  If  x*  is  optimal  in  (P)  and  any  of  the  following  hold, 

(a)  (interiority)  xo  €  n(.Yo)  (ri  is  relative  interior)  and  intf-YjOFf )  ^  0  relative  to  aff(.Y,U>', ): 
or, 

(b)  (slack  time)  for  any  t,  xt  €  XtC\Yt  there  exists  {xt+i,...},  Ft(xt)  =  inf  /r(xr,  xr+j ).  T; 
t  such  that  the  common  resource  groups  required  for  any  set  of  parts  are  idle  at  T:  or 

(c)  (asymptotic  penalty  free  schedule)  for  any  t,xt  6  -Yjny),  there  exists  T  >  <,{xj.x}+1....} 
such  that  /i(i|,xj+1)  +  ELi/^r^r+i)  <  oo,  /r(x}.,xj.+i)  =  0  for  all  r  >  T\  or 

(d)  (decreasing  fixed  matchup  cost)  there  exists  a  feasible  solution  to  P,  x'  =  {x(,xj+1 _ } 

such  that  for  any  t,x<  6  X t  n  Yt,Tx,K  =  1,2,. . .  ,Tk  <  Tk+ i,  such  that 

tk-\ 

^2  /(xr,x,+j) +  /T|f(XTJf,x'7W+1)  <  oo 

and  limft-^oc  frK (xtk,x'tk+\)  =  0; 
then  there  exists  p*,  t  =  0,1,...  such  that 

(i)  F\x;)  -  p\x;  <  Fi(xt)  -  p(*xf  for  all  x,  €  X„t  =1,2,... 

(ii)  /i(x;,x('+1)  -  PtXf  +  Pt+iXt+i  <  /«(x«,xt+j) -PtXt  +p;+1x,+i,  for  all  (xt,x«+i)  such  that 
/i(x,,x,+i)  <  oo. 

Proof:  See  Appendix. 

The  asymptotic  stability  of  solutions  follows  again  from  the  structure  of  the  objective  function. 

The  following  theorem  states  that  optimal  solutions  from  varying  initial  conditions  asymptotically 
approach  each  other.  That  is,  the  eventually  optimal  path  is  insensitive  to  the  initial  conditions. 


Let  Z *  be  the  set  of  solutions  (rt,ct+i)  such  that 

ftixt  i x  t+\)  ~~  Ptxt  "b  Pi+l^t+l  =  ftizti  2«+l )  —  Pt  st  "b  Pt+\‘-t+ 1  • 

Our  result  is  that  all  optimal  solutions  approach  Z*  regardless  of  initial  condition. 

Theorem  2:  Let  x*  be  optimal  for  (P)  with  initial  condition  xo,  price  support  p*  and  facets  Z* 
defined  as  above.  Let  x'  be  optimal  in  ( P )  with  initial  condition  x'0  and  price  supports  p'.  Then, 
for  any  e  >  0,  there  exists  T  <  o o,  such  that,  for  all  t  >T, 


inf 


(x4,xt+j)  (zt,  Zt+l , 


<  €. 


(2) 


Proof:  See  Appendix 

Beyond  the  asymptotic  approach  to  a  turnpike  solution,  we  can  also  show  that  matchup  to 
a  turnpike  occurs  as  quickly  as  possible  if  the  only  objective  is  to  minimize  tardiness  as  defined 
above.  This  result  further  justifies  our  use  of  the  matchup  approach. 


Theorem  3:  Assume  in  Problem  (P)  that  =  0  for  all  i  and  that  x*  is  an  optimal  solution 
given  xjj.  Let  x'0  <  xj  be  an  alternative  initial  state.  If  there  exists  a  feasible  solution  x  such  that 
xo  =  Xq  and  x(  =  x*  for  some  t  <  oo,  then  there  exists  an  optimal  solution  x'  with  initial  condition 
Xq  such  that  x\  =  x*. 

Proof:  See  Appendix. 

The  impact  of  this  result  is  that  the  optimal  schedule  beginning  at  the  disruption  state  matches 
up  with  an  optimal  preschedule  as  soon  as  possible  with  finite  cost.  Two  questions  arise:  how  fast 
can  this  matching  up  take  place,  and  what  happens  if  additional  disruptions  occur  in  the  system? 

The  answer  to  the  first  is  dependent  on  system  parameters.  For  example,  a  system  with  well 
distributed  ample  slack  resources  can  matchup  quickly.  Note  that  this  observation  leads  to  the 
important  question  of  how  the  system  and  preschedule  can  be  designed  to  make  rescheduling  easier. 

The  second  question  results  in  an  assumption  that  disruptions  are  spaced  far  enough  apart 
that  match  up  is  possible  between  disruptions.  This  situation  occurs  in  the  real  test  problems 
studied  below.  If  this  assumption  is  not  valid,  the  system  is  so  disrupted  that  no  preschedule  will 
be  useful  and  full  dynamic  scheduling  is  likely  near  optimal. 


As  an  implementation  issue,  we  wish  to  use  this  matchup  philosophy  to  reschedule  a  system 
upon  disruption.  We  need  to  determine  the  matchup  time  and  the  optimal  schedule  from  the 


disruption  time  to  the  matchup  time.  In  practice,  neither  is  a  simple  task.  Below  we  describe  a 
heuristic  born  of  the  matchup  philosophy  and  test  it  against  real  problems  from  the  automotive 
industry. 

We  define  a  job  as  production  unit  (e.g.,  part)  with  a  single  due  date  and  release  date  that 
cannot  be  subdivided  for  processing.  A  lot  is  a  collection  of  jobs  from  the  same  part  type  that 
does  not  require  setups  between  jobs.  A  tool  is  a  resource  required  in  addition  to  a  machine  for 
processing  a  job. 

Assume  that  a  preschedule  has  been  constructed  and  implemented.  When  a  disruption  occurs 
we  seek  to  reschedule  jobs  on  machines  and  tools  to  minimize  total  weighted  tardiness  such  that 
no  release  dates  are  violated. 

Since  this  problem  is  NP  hard  (even  rescheduling  a  single  machine  subject  to  release  dates  to 
minimize  total  tardiness  is  unary  NP  hard,  Graham,  et.  al.  [1079]),  we  use  a  heuristic  algorithm. 
We  first  attempt  to  match  up  on  individual  machines.  If  this  is  unsuccessful  we  reallocate  lots 
across  machines  and  resequence  jobs  on  individual  machines.  This  allocation  and  sequencing  step 
is  similar  to  that  in  Dogramaci  and  Surkis.  In  the  algorithm  given  below  note  that  Tl,TmM,DT 
and  EPS  are  user  defined  parameters. 

Matchup  Scheduling  Algorithm  (MUSA) 

Step  0  :  For  each  disrupted  machine,  set  a  minimum  match-up  time,  IT.  Let  T  =  IT.  Go  to  Step  1. 

Step  1  :  On  each  disrupted  machine,  resequence  all  jobs  scheduled  before  T.  Evaluate  the  schedule 

COST  on  each  disrupted  machine.  If  COST  <  EPS  on  all  machines,  STOP.  Else,  go  to  Step 

2. 

Step  2  :  Let  T  =  T  +  DT.  If  T  >  TmhX,  go  to  Step  3.  Else,  go  to  Step  1. 

Step  3  :  Expand  the  set  of  machines  to  be  rescheduled  to  include  all  machines  that  share  job  com¬ 

patibilities  with  the  current  set  of  disrupted  machines.  Reallocate  lots  across  these  machines. 

Go  to  Step  0. 

The  procedures  used  for  the  allocation  and  sequencing  steps  are  described  below  and  include 
the  allocation  of  a  finite  set  of  tools.  The  following  sections  describe  the  alternatives  implemented 
for  these  procedures. 


3.  Single  Machine  Sequencing 


In  Step  1,  MUSA  heuristically  reschedules  jobs  on  a  single  disrupted  machine.  This  method 
must  run  quickly  due  to  the  real  time  environment.  The  goal  of  the  procedure  is  to  minimize 
the  total  cost  for  late  processing  of  the  jobs  where  each  job  is  started  no  earlier  than  its  release 
date.  To  increase  scheduling  flexibility,  lots  are  broken  into  jobs.  The  overall  objective  is  to  find  a 
sequence  of  start  times  to  minimize  the  sum  of  setup  and  tardiness  costs  without  violating  release 
date  or  tool  availability  constraints. 

The  heuristic  calculates  six  feasible  schedules  based  on  different  ordering  rules  and  chooses 
that  with  least  cost.  The  ordering  rules  are  1)  a  shortest  processing  time  ordering  (SPT),  2)  an 
earliest  due  date  ordering  (EDD),  3)  a  modified  due  date  ordering  (MDD),  4)  a  priority  index 
ordering  (API),  5)  a  ratio  rule  and  6)  the  ordering  based  upon  the  given  sequence.  The  MDD 
rule  was  taken  from  the  heuristic  developed  by  Baker  and  Bertrand  [1982].  The  API  rule  was 
taken  from  the  heuristic  developed  by  Morton  and  Rachamadugu.  The  ratio  rule  is  based  on  a 
comparison  of  the  remaining  processing  time  of  a  job  to  the  length  of  time  available  until  the  job  is 
due.  The  given  sequence  ordering  rule  yields  a  pushback  of  the  schedule  in  use  when  a  disruption 
occurs. 

For  details  on  the  implementation  of  the  single  machine  sequencing  heuristic  see  Bean.  Birge. 
Mittenthal  and  Noon  [1986], 

4.  Multimachine  Lot  Reassignment 

If  the  best  single  machine  solution  results  in  excessive  overtime  or  premium  freight  charges, 
the  multimachine  lot  reassignment  program  is  used  to  redistribute  lot  to  machine  assignments. 
The  focus  of  the  reassignment  problem  is  to  determine  a  feasible  lot  schedule  across  several  ma¬ 
chines.  The  multimachine  reassignment  uses  lots  for  rescheduling  in  order  to  keep  the  problem 
size  manageable  and  reduce  additional  setups.  We  examined  two  approaches  for  reassignment:  a 
multiple  choice  integer  program  (MCIP)  formulation  solved  using  the  technique  of  Bean  [19S4i. 
and  a  priority  rule  dynamic  assignment  heuristic. 

The  reassignment  procedure  can  shift  lots  across  machines  to  create  a  feasible  multimachine 
schedule  without  additional  setups.  The  single  machine  resequencing  algorithm  can  then  be  applied 
to  order  jobs  to  achieve  further  penalty  reductions. 

4.1  Integer  Programming  Approach 

The  multiple  choice,  zero  one  program  has  decision  variables  xtJ,  which  equal  1  if  lot  i  is 


assigned  to  machine  j  and  equal  0  otherwise.  If  lot  i  can  be  run  on  some  subset  of  the  machines. 
Mi,  with  cardinality  n,  then  n  variables  are  created  for  lot  i.  To  ensure  that  lot  i  is  scheduled  on 
exactly  one  machine,  the  logical  constraint 

£  *ij  =  1  (3) 

j£Mi 

is  used.  The  model  also  adds  machine  utilization  constraints  to  ensure  that  the  scheduled  processing 
time  does  not  exceed  the  available  processing  time.  These  constraints  are  written 

^PijXij  <Hj,  (4) 

where  pi}  is  the  processing  time  of  lot  i  on  machine  j  and  H}  is  the  available  processing  time  on 
machine  j. 

Additional  constraints  for  feasibility  of  the  lot  to  machine  assignments  are  needed  to  prevent 
the  assignment  of  more  than  one  lot  to  one  machine  simultaneously  and  to  prevent  the  use  of 
a  single  tool  on  two  machines  simultaneously.  A  lot’s  window  is  the  interval  between  its  release 
date  and  due  date.  Ideally,  we  would  like  to  allow  any  placement  within  this  window.  For  the 
integer  programming  approach,  however,  the  resulting  formulation  is  too  complicated  to  solve  in 
real  time.  To  simplify  this  formulation  placements  within  the  window  were  discretized.  One  or 
more  possible  placements  are  determined  for  each  lot  to  machine  combination.  The  formulation 
is  then  altered  so  that  x,j  =  1  if  lot  i  is  placed  in  position  j  (which  has  an  associated  machine). 
This  increases  the  number  of  variables,  but  the  structure  of  the  MCIP  is  retained  allowing  efficient 
solution.  One  approach  allows  three  placements,  justified  left,  centered,  or  justified  right.  The 
other  fixes  placements  to  avoid  conflicts.  For  details  see  Bean,  Birge,  Mittenthal  and  Noon. 

Given  start  times,  Si,  and  finish  times,  Fi,  for  all  potential  lot  placements,  constraints  are 
created  to  avoid  the  assignment  of  substantially  overlapping  lots  to  the  same  machine.  Some  overlap 
is  allowable  as  it  can  be  resolved  later  by  single  machine  sequencing  or  overtime.  The  allowable 
overlap  is  denoted,  RELAX.  This  relaxation  helps  compensate  for  modeling  the  continuous  lot 
start/finish  by  discrete  assignments.  Constraints  are  added  such  that  on  machine  k.  if  F,  > 
(Sj  +  RELAX)  then 

X ik  4"  Xjh  ^  1,  ^  ■ 

The  final  set  of  constraints  prevents  the  simultaneous  scheduling  of  two  lots  which  use  the 
same  tool.  If  lot  i  on  machine  j  and  lot  k  on  machine  /  use  the  same  tool  and  F,  >  5*  and  F*  >  >, 


then 


Xij+Xkt<l.  (6) 

Constraints  (3),  (4),  (5),  and  (6)  provide  the  basis  for  the  multimachine  reassignment  program. 
With  this  foundation,  the  problem’s  objective  can  take  on  a  number  of  forms.  We  tested  two 
alternatives:  minimize  the  number  of  lot  to  machine  assignment  changes,  and  maximize  the  sum 
of  squared  processing  times  scheduled.  A  weighted  combination  of  these  objectives  produced  the 
best  results. 

Other  considerations  incorporated  in  the  program  include  constraints  on  machine  utilization  to 
balance  machine  workload  and  variable  machine  production  rates.  Other  possible  additions  include 
constraints  on  manpower  requirements  for  discrete  intervals  to  inhibit  manpower  imbalance,  and 
adherence  to  precedence  requirements. 

4.2  Priority  Rule  Approach 

An  alternative  approach  to  reassigning  lots  to  machines  via  the  MCIP  solution  is  to  use  a 
priority  rule  strategy.  The  alternative  we  describe  is  similar  to  the  approach  in  Dogramaci  and 
Surkis  and  to  that  of  the  single  machine  resequencing  heuristic.  Both  develop  initial  schedules 
based  on  a  number  of  different  rules  and  then  select  one.  The  heuristic  develops  an  initial  lot 
sequence  over  the  machine  group  and  submits  them  one  at  a  time  for  resequencing  to  the  single 
machine  heuristic.  A  description  of  this  heuristic  algorithm  is  given  ;n  Bean,  Birge,  Mittenthal, 
and  Noon. 


5.  Experimental  Results 


The  matchup  scheduling  algorithm  (MUSA)  including  the  alternatives  for  lot  allocation  and 
job  sequencing  was  coded  in  FORTRAN  and  implemented  on  an  IBM  4381  computer.  The  program 
was  applied  to  a  set  of  problems  from  an  automotive  manufacturer.  The  goal  of  the  comparisons 
on  this  data  set  was  to  determine  the  value  of  the  matchup  approach  relative  to  the  simple  static 
and  dynamic  approaches,  to  evaluate  the  use  of  the  MCIP  allocation  scheme  relative  to  priority 
rules  and  to  determine  the  relative  values  of  the  individual  sequencing  rules. 


The  problem  set  consisted  of  eight  disruption  scenarios  from  a  facility  with  two  machines,  and 
five  disruption  scenarios  from  a  facility  with  ten  machines.  Data  for  Facility  1  included  58  lots  (250 
jobs)  scheduled  over  two  fully  compatible  machines  with  an  average  utilization  of  76 7c.  Facility  2 
consisted  of  25  lots  (203  jobs)  scheduled  over  ten  partially  compatible  machines  with  an  average 


utilization  of  46%.  The  disruptions  were  chosen  to  represent  common  difficulties  which  render  a 
preschedule  infeasible:  machine  breakdowns,  tool  unavailabilities,  release  or  due  date  changes,  and 
order  quantity  increases.  The  preschedules  were  developed  by  the  manufacturer  and  contained 
some  unresolvable  tardiness  before  the  disruptions  were  added.  This  was  due  to  earlier  disruptions 
that  forced  ready  times  plus  processing  times  to  be  greater  than  due  dates.  This  inherent  tardiness 
was  a  lower  bound  on  the  total  tardiness.  Tables  la  and  lb  identify  the  type  of  disruption  and 
gives  the  inherent  tardiness  lower  bounds. 

Parameters  were  chosen  to  reflect  operating  conditions  in  the  facilities.  The  solution  horizons 
generally  included  about  70  percent  of  the  lots  and  were  chosen  so  that  the  cumulative  idle  time 
was  at  least  twice  the  average  lot  processing  time.  The  MC1P  formulation  initially  allows  no 
overlap  ( RELAX  =  0)  among  lot/machine  assignments.  If  no  feasible  integer  solution  can  be 
found,  RELAX  is  successively  increased  until  enough  constraints  are  relaxed  that  a  solution  can 
be  found. 

Four  strategies  for  rescheduling  were  evaluated  to  illustrate  the  nature  of  the  matchup  prob¬ 
lem.  Strategy  1  gives  a  static  solution.  When  a  disruption  occurs,  the  machine  assignments  and 
job  sequences  stay  the  same,  only  the  job  start  and  finish  times  are  shifted  to  accommodate  the 
disruption.  Strategy  2  gives  a  fully  dynamic  myopic  solution  using  the  myopic  priority  rule  heuris¬ 
tic.  For  this  strategy  the  reassignment  heuristic  wa  used  to  reschedule  lots  according  to  one  of 
three  selection  rules:  EDD,  MDD  and  LWS.  The  last  rule  orders  lots  based  on  least  slack  time  in 
its  window. 

Strategy  3  considers  only  partial  look  ahead.  For  this  strategy,  the  multimachine  reassignment 
problem  was  modeled  as  an  MCIP.  The  preschedule  lot/machine  assignments  and  starting  times, 
however,  were  excluded  from  the  formulation.  Rather  than  dynamically  scheduling  forward  through 
time,  as  in  the  priority  rule  heuristic,  this  approach  solves  the  problem  over  the  matchup  horizon 
and  provides  some  degree  of  “look  ahead”  when  considering  the  impact  of  scheduling  a  lot.  After 
the  MCIP  determines  the  multimachine  lot  schedule,  the  single  machine  resequencer  is  used  to 
obtain  further  penalty  reductions. 

Strategy  4  represents  full  matchup  rescheduling  using  the  MCIP  with  the  preschedule  job 
assignments  included  in  the  formulation.  This  strategy  displays  the  sirength  of  combining  the  look 
ahead  aspects  of  the  MCIP  with  the  option  of  returning  to  the  preschedule  when  it  is  advantageous 
to  do  so.  After  the  MCIP  determines  the  multimachine  lot  schedule,  the  single  machine  resequencer 
is  used  as  in  the  previous  strategy. 


The  tardiness  results  for  the  four  strategies  are  displayed  in  Tables  la  and  lb.  The  values 
represent  total  job  tardy  hours  across  all  machines.  The  tardiness  results  for  Strategy  1  illustrate 
the  penalties  incurred  when  machine  compatibilities  are  not  utilized  during  rescheduling.  When  the 
preschedule  is  pushed  back  to  accommodate  the  disruption,  only  the  jobs  on  the  disrupted  machine 
are  affected.  Although  this  strategy  preserves  the  preschedule  sequence  and  machine  assignments, 
it  is  limited  since  jobs  may  not  be  offloaded  to  compatible  machines. 

In  Strategy  2,  among  the  three  selection  rules  tested,  EDD,  MDD  and  LWS,  the  LWS  rule 
had  the  lowest  average  tardiness,  however,  it  did  not  significantly  outperform  the  others.  The 
tardiness  comparison  between  Strategy  2  and  Strategy  3  is  mixed.  The  partial  look  ahead  approach 
performed  better  on  the  Facility  1  problems  but  markedly  worse  on  the  Facility  2  problems.  This  is 
due  to  MCIP’s  difficulty,  without  the  preschedule  assignments,  in  finding  a  feasible  schedule  across 
the  many  machines  of  Facility  2. 

Strategy  4  yielded  the  least  tardiness  for  both  facilities  among  all  strategies.  By  including  the 
preschedule  assignments  in  the  MCIP  formulation,  the  model  switched  job/machine  assignments 
only  as  needed  to  correct  for  the  disruption.  The  preschedule  assignments  then  served  as  a  good 
completion  to  the  solution  schedule. 

The  appropriateness  of  a  rescheduling  strategy  must  also  consider  machine  assignment  changes 
and  computation  times.  Table  2  displays  the  number  of  job/machine  assignment  changes  for 
Strategies  2,  3  and  4.  Since  the  simple  preschedule  pushback  approach  of  Strategy  1  does  not 
move  jobs  across  machines,  no  values  are  given.  For  the  Facility  1  problems,  Strategy  2  yielded  a 
significantly  higher  number  of  job  changes  than  Strategies  3  or  4.  This  is  to  be  expected  since  the 
reassignment  heuristic  does  not  attempt  to  maintain  original  assignments  as  the  MCIP  does.  All 
three  strategies  had  few  job/machine  changes  for  the  facility  2  data  since  many  jobs  had  only  one 
compatible  machine. 

Table  2  also  displays  the  computation  times  for  Strategies  2,  3  and  4.  The  time  for  Strategy 
1  was  negligible  and,  hence,  omitted.  The  Strategy  2  heuristic  has  the  advantage  of  being  very 
fast  when  compared  to  the  MCIP  run  times  of  Strategies  3  and  4.  The  latter  strategies  require 
the  solution  of  MCIP’s  with  up  to  225  variables  and  several  hundred  constraints.  As  the  number 
of  lots  increases,  the  CPU  times  for  strategies  1  and  2  can  be  expected  to  increase  approximately 
linearly  while  those  of  strategies  3  and  4  can  be  expected  to  increase  more  sharply. 

Table  2  includes  the  final  RELAX  values  for  Strategies  3  and  4.  This  represents  the  amount 
of  job  overlap  that  had  to  be  allowed  before  a  feasible  integer  solution  could  be  found  by  the  M<  IP. 


The  Strategy  3  RELAX  averages  are  higher  than  those  of  Strategy  4  indicating  greater  difficulty 
in  finding  a  feasible  solution.  The  direction  from  prescheduling  assignments  not  only  resulted  in 
lower  RELAX  values,  but  also  lower  average  run  time. 

The  next  set  of  comparisons  evaluates  the  matchup  scheduling  using  the  MCIP  solution  relative 
to  the  priority  rule  reassignment.  The  priority  rule  heuristic  was  run  with  each  of  its  five  rules. 
The  resulting  schedules  were  then  evaluated  with  respect  to  lot  tardiness  and  the  minimum  tardy 
schedule  was  chosen.  The  single  machine  resequencer  was  applied  to  the  chosen  schedule  to  evaluate 
and  improve  the  job  tardiness.  This  approach  was  compared  to  Strategy  4  which  used  the  MCIP 
for  lot  reassignment.  Table  3  displays  the  results  of  the  comparison. 

For  the  heuristic  selection  rule  choice,  the  results  show  a  fairly  even  distribution  among  the 
five  rules  tested  and,  on  several  problems,  identical  lot  schedules  among  some  of  the  rules.  In 
terms  of  mean  tardiness,  the  MCIP  approach  displays  a  slight  edge  over  the  heuristic.  This 
difference  arises  from  several  problems  in  which  the  heuristic  performed  quite  poorly  compared  to 
the  MCIP  approach  because  of  the  myopic  versus  look  ahead  characteristics  of  the  two  approaches. 
The  average  number  of  machine  changes  also  favors  the  MCIP,  especially  among  the  Facility  1 
problems.  The  heuristic,  however,  displays  much  lower  CPU  times  than  the  MCIP. 

6.  Conclusions 

This  paper  presents  a  framework  for  scheduling  production  facilities  when  disruptions  inval¬ 
idate  preplanned  schedules.  It  is  shown  that  if  disruptions  are  sufficiently  spaced  over  time,  that 
the  optimal  rescheduling  strategy  is  to  match  up  with  the  preschedule.  Heuristics  are  designed 
to  implement  this  philosophy  for  problems  with  multiple  resources,  setups  and  release  dates.  The 
methods  were  tested  on  a  practical  set  of  test  problems  from  an  automobile  manufacturer.  The 
results  showed  that  the  matchup  scheduling  approach  provided  significantly  better  results  than 
the  pure  static  and  dynamic  strategies.  The  experiments  also  indicated  that  the  MCIP  integer 
programming  solution  to  the  allocation  problem  may  be  useful  when  utilizations  are  high. 
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APPENDIX 


Theorem  1:  If  x*  is  optimal  in  (P)  and  any  of  the  following  hold, 

(a)  (interiority)  x0  €  n(.Y0)  (rt  is  relative  interior)  and  int(.Ytny<)  ^  (J  relative  to  aff(.Y<U}'< ): 
or, 

(b)  (slack  time)  for  any  t,xt  G  -YtnF(  there  exists  {xt+J,...},  Ft(xt)  =  inf  fr(xr- Xr+\ ;.  T 
t  such  that  the  common  resource  groups  required  for  any  set  of  parts  are  idle  at  T ;  or 

(c)  (asymptotic  penalty  free  schedule)  for  any  t,xt  G  -Yt  H  Y<,  there  exists  T  >  t,  [x),x]+1 _ } 

such  that  /«(xt,xj+i)  +  YJT=t  fr(x\.,x *+,)  <  oo,  /r(x).,xJ.+ j)  =  0  for  all  r  >  T\  or 

(d)  (decreasing  fixed  matchup  cost)  there  exists  a  feasible  solution  to  P,  x'  =  {xi,x'l+1 _ } 

such  that  for  any  t,xt  €  -Yt  n  \\,Tk,K  =  1,2 , ...  ,7V  <  Tk+  i,  such  that 

Tk- i 

/(^r,Xr+l)  +  fTK(XTK,XTK  +  i)  <  00 

and  hmA-^x/^fxr^XT^j)  =  0; 

then  there  exists  p*,t  =  0, 1, . . .  such  that 

(i)  Ff(x*)  -  p‘txmt  <  Pt(xt)  -  PtXt  for  all  xt  €  Xt,t  =  1,2,... 

(ii)  ft(X(,xmt+1)  -p?xj  +p?+iX*+1  <  /«(x,,xt+1) -p*x,  +p*+iXl+,,  for  all  (x,,x,+1 )  such  that 
ft{Xt,  Xf+] )  <  oo. 

Proof:  McKenzie  [Lemma  1,  1976]  proves  the  result  given  (a).  This  condition  may  not  hold, 
however,  in  the  scheduling  context  for  x0  €  dXo  and  the  interiority  of  .Y,  n  Yt  may,  indeed,  be 
difficult  to  verify.  The  other  conditions  are  reasonable  assumptions  that  may  be  more  readily 
verified.  We  show  that  each  implies  (i)  and  use  (i)  and  the  structure  of  ft  to  show  (ii). 

Condition  (b)  implies  that  there  is  sufficient  slack  in  the  schedule  for  all  resources  for  common 
parts  to  be  free  eventually.  This  condition  should  hold  for  all  but  heavily  loaded  systems.  Note 
that  the  objective  in  (P)  is  completely  separable  among  sets  of  parts  sharing  common  resources  so 
we  can  assume  without  loss  of  generality  that  all  parts  have  a  common  resource  and  all  resources 
are  eventually  simultaneously  idle.  We  wish  to  show  that  F*(  c<)  is  subdifferentiable  at  x(. 

Consider  Xj  G  X(  n  Yt,Fi{xt)  =  fr(xr,  xr+\ )  and  x'<  =  xt  +  6.  The  idleness  property 

implies  that  there  exists  a  feasible  schedule,  x",  such  that  x"  =  x't  and  x'r'  =  xr,r  >  T.  We 
claim  there  exists  an  optimal  schedule  from  x't  among  x",  r  >  T ,  i.e.,  there  exists  {x'<+ ,....}  such 


m 


that  Ft(x't)  =  E^t/r(I'r>xr+i)’xr  =  xr,  r  >  T,  for  all  ||<!>||  <  €  and  some  e  >  0.  Suppose 
not,  then  there  exists  some  processing  after  T  in  x<  that  can  be  moved  before  T  to  reduce  cost 
or  processing  before  T  in  xt  that  can  be  moved  after  T  to  reduce  cost.  In  either  case,  the  same 
change  in  processing  can  be  moved  to  T  with  reduced  cost  by  the  structure  of  This  implies 
that  {xr,r  >  f}  is  not  optimal,  a  contradiction. 

From  the  above,  only  a  finite  number  of  states  xr  are  changed  in  an  optimal  path  from  x,  to  an 
optimal  path  from  x't.  Since  ft  has  bounded  slope,  there  exists  A'  <  (T-t)  £"=1  max{u’,,  u, }  <  x 
such  that  Ft(x't)  >  Ft(xt)  -  A'||xJ  —  x(j|.  Hence,  Ft  is  subdifferentiable  at  x(,  proving  (i). 

A  similar  argument  is  used  if  (c)  holds  by  noting  that  only  a  finite  number  of  costs  are  reduced 
in  an  optimal  path  from  x't  from  an  optimal  path  from  xt.  This  again  implies  subdifferentiability. 

Condition  (d)  can  be  interpreted  as  a  generalization  of  (c),  in  which,  a  finite  cost  path  is 
eventually  ootainable  from  every  feasible  path  at  decreasing  cost.  Note  that  /Ur.Xr-H )  + 

/tk(xTk  ^tk+i)  h*3  a  finite  number  of  terms  with  bounded  slope  and  is  hence  subdifferentiable. 
Hence,  there  exists  pi'  such  that 

Tk-1 

52  f(*r,x r+i)  +  /T*(xrk<xr*+i)  -P^t 

r=l 

7* -1 

<  52  +  frK  ( 1 7V  ’  xTk  + 1  )  -  Pt'x ,  I  ~ 


for  all  x"  €  Xt  .  Rewrite  (7)  as 


Pt(xt  xl)  ^  52  /(xrixr+l)  ~  52  /(Ir'  xT+\ )  + 


fTK(xTK,xTK  + 1)  ~  fTK(xTK,xTK  + 1)- 


Note  that  |/r‘(x()|  <  oo  and  <  oo.  Hence  pi'  has  a  limit  point,  p<,  and  by  (d; 


Pt(xt~x")  <F‘(x'r')  -F‘(x«). 


for  all  x'r'  €  A(. 
Note  that 


F^Wi)  =  inf  {/.-dx.-M;)  +  F*(x'«)} 
-  /*-! (x,_i ,  x()  +  F*(x»). 


V w 


TO' 


Consider  the  function,  g(x't)  defined  by 


g(A)  = 

-  fi-l(x't-l'X't)  +P?-l*i-l  (11) 

<  F\x't), 

for  any  (xj_,,i',)  such  that  /<-i(x't_,, xj )  <  c©  and  note  that  g  is  also  subdifferentiable  for  x't  6  Yt. 
Hence,  there  exists  p't  such  that 

g(x't)  -  p',x't  <  g(x't)  -  ptx"  <  F‘(x")  -  p'tx't',  ( 12) 

for  all  x"  £  Xt  n  Yt.  Now,  let  x','  =  x*,  and  p\  =  p4  to  obtain 

ft—  l(X(_j,Xt)  —  P(_jX,_j  —  ft—  l(xt_  Xj)  +  Pf_iX,_j  —  p,  x,  ^  — p,  X,  ,  (  1  3 1 

for  all  (xj_i , x't)  feasible,  proving  (ii).« 


Theorem  2:  Let  x4  be  optimal  for  ( P )  with  initial  condition  xo,  price  support  p*  and  facets  Z* 
defined  as  above.  Let  x'  be  optimal  in  ( P )  with  initial  condition  x’0  and  price  supports  p'.  Then, 
for  any  e  >  0,  there  exists  T  <  oo,  such  that,  for  all  t  >  T, 


.  ‘of  -  (*«,*«♦!  )||  <  «• 

Proof:  Let  x'  have  supporting  prices  p'.  From  (ii)  of  Theorem  1,  we  have  for  any  x*  £  Z4, 
(x<  —  X|  )  —  P*+i(x(+,  —  xf+1 }  <  /|(X,,X,+  1)  —  /«(x,,  X,+  , ) 

—  Pt(xt  -  zt)  ~  Pt+l(xl+l  —  ‘(+1  )• 

Let  t’c(x,4)  =  {pt  -  p't)(x't  -  x4)  .  Inequality  (15)  implies  that 


ill) 


(15) 


r,(x;+l)<r,+,(x,4}  (lfi) 

for  all  t  and  any  (x,*,x4+1)  £  Z*.  By  summing  inequalities  (ii),  for  all  7\ 

T 

Pt+i  ( XT+ 1  —  XT+ 1 )  —  ^  '(/e(xt  >  x(+i )  —  /i(Xj,  x,+j ) )  —  Pq(  Xo  —  x0 ) ,  ( ,  t  : 

oo 


and 


T 

Pr+1  (xT+1  _  i7'+i)  —  ^^(/l(x(.  It+ 1  )  -  ft(x(  ,  xl+ j ))  -  p0(x0  -  Xo  ). 
oo 


c 
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From  the  finiteness  of  F°(x')and  F0(x*),  both  right  hand  sides  in  (17)  and  (18)  are  uniformly 
bounded  for  all  T.  Hence,  we  can  assume 


t’f ( £j  )  >  K  >  -00, 

for  all  (  and  x*. 

Now,  if  (14)  does  not  hold,  then,  for  all  T,  there  exists  some  T'  >  T  such  that 

.  inf.c7.  ll*'«  - --<11  >  <• 


If)  | 


(20) 


(21 


By  boundedness,  let  the  infinium  in  (18)  be  attained  at  z*.  Then,  by  the  structure  of  1  =  V . 

A(:;,s;+i)-p;*;+p;+i*;+i 

—  /<(x,,  xt+I )  ~  Ptxt  +  p<+ix<+i  —  7 IK :t  i  -<+i )  ~  (xf  xt+i ) I! 

where  ■>  >  min(ut,,u,)  >  0.  By  definition  of  Z*,  for  any  (x*,x*+1) 

fi(xt  >x<+i )  ~  Pt  xt  +  Pt+ixt+i  —  ft(:t '  rt+t )  ~  Pt:t  +  P<+i  r<+i  • 

From  (21)  and  (22),  we  have  for  any  T, 

T  T 

y!  /»(x<  -  x<+i )  ~  ^  /<(x<-  x<+i ) 


l»0 


1-0 


<  (PS  -  Po)(xo  -  xo)  -  (PT  -  Pt)(x’t  -  xr)  -  £  7ll(x«* -,*+l)  -  « 


J! 


r«o 


By  ( If)),  if  ( 14)  does  not  hold,  then  the  right-hand  side  of  (23)  approaches  -oc  as  T  approaches 
so.  This  contradicts  the  finiteness  of  /^(x')* 


Theorem  3:  Assume  in  Problem  ( P )  that  u,  =  0  for  all  i  and  that  x*  is  an  optimal  solution 
given  Xq.  Let  x'0  <  xj|  be  an  alternative  initial  state.  If  there  exists  a  feasible  solution  x  Mich  that 
x0  =  Xq  and  xt  =  x*  for  some  t  <  oc,  then  there  exists  an  optimal  solution  x'  with  initial  condition 
Xq  such  that  x'(  =  xj. 

Proof:  First  construct  a  feasible  solution  x  such  that 

(i)  Jr(t)  =  x'r(0  if  x'r(t)  <  x*(«), 

(ii)  xr(j  )-  xr_|  (j  )>  x*(; )-  x*_,(j )  for  all  j  sharing  resources  with  i  if  x'ri  i  >  r!  -  and 


A  solution  satisfying  (i)  and  (ii)  satisfies  xr  <  x*.r  =  1,2 . f,  since  no  processing  is  made 

beyond  x*(i)  for  each  i.  The  feasibility  of  a  path  from  x'0  to  x*  only  requires  5Z"_,(xo(0  -  x()((;; 
available  capacity  units  on  the  resources  required  by  those  parts.  This  processing  can  be  completed 
anywhere  in  [0 , /J  such  that  xo  >  Xq  implies  r(j)  <  0  for  all  i  such  that  xo(?)  <  xj, ( * )  and  no  other 
resource  constraints  are  present.  Therefore,  x  satisfying  (i)  and  (ii)  can  be  extended  to  satisfy 
X|  =  x*  and  xr  =  x* ,  r  =  0, 1 , . . . ,  /. 

Suppose  x  ^  x'  for  any  x'  optimal  given  xj,.  Then,  there  exists  a  sequence  Ax  =  (Axn  = 
0,  Axi ,  Axj, . . .)  such  that  x'  =  x  +  Ax  and  F(x  4-  Ax)  <r  F(x).  We  will  show  that  this  implies 
that  there  exists  a  sequence  Ax*  such  that  F(x*  +  Ax*)  <  F(x*).  To  show  this,  let  F(x')  = 
F,(x'(«)),  where  F,  is  the  contribution  to  the  overall  objective  value  from  part  i.  Suppose 
F,(x'(i))  <  F,(x(j)).  Then  there  exists  r  such  that  Axr(i)  -  Axr_i(t)  >  0.  For  each  such  r.  note 
that  x,(i)  =  x*(i),  so  x*(i)  +  Axr(i)  =  xr(i)  +  Axr(r)  and 

F,(x(i)  +  Ax(t))  -  F,(i(i))  =  F,(x*(i)  +  Ax’(i))  -  F,(x*(0).  (-4; 

where  Ax*(i)  =  Ax(i). 

In  order  for  x(i)  +  Ax(i)  to  be  feasible,  however,  there  must  exist  Ax,(j)  -  Axr-i(;'  ‘I  *1  for 
parts  j  requiring  common  resources  with  i.  For  each  such  j,  we  can  define  Ax*  recursively  from 
r  =  1  bv 

AxJO)  -  Ax;_,0)  =  min {x'(j)  -  x*_j(j),Ax*(j)  -  ArJ., (;)}.  (25. 

By  (ii),  this  ensures  the  feasibility  of  x*  +  Ax(i)  +  £  6y(l)  Ax*(; ),  where  J(i)  =  (set  of  parts 
sharing  resources  with  t}.  For  each  j  €  ./( r ) ,  note  that 

F,(x(;))  ~  F(x(;)  +  Ax(j))  =  Wj^Xr(j) 

Jr|t,(j)<0) 

■f  ^  U’j(xr(j)  +  Axr(j))  +  ^  n.  ijn 

{rltr(j)>0,tr(j)  +  tlxr(j)<0)  f  rjt ,  (  J  )>0 .1 ,  (  J  )  +  r  <  J  )>0  } 

Since  xr  <  x*  and  Ax*  <  Axr,  from  (26), 

Fj(x(j  ) )  -  Fj(x(j)  +  Ax(  j  i  i  <  Fjix'i  j))  -  Fjlx’i  j :  f  A x*(;)).  \'2~ 


Fonstruct  Ax*  using  this  procedure  for  each  i  such  that  F,(x'(i))  <  F,(x(i!).  It  follow-  from 


i, 


w 


*  £  Fi(x7‘)  +  Ax*(i))-^(**(e))+  ]T  fj(x*0')  +  Ax*(»j 

/={«|Fi(x'(i))<F.(i(»))}  {j€  •),>€/) 

<^F,(x'0'))-irI(x(0)+  £  ^(x'O))  -  Fj(x(j))  <  0, 

but  this  contradicts  the  optimality  of  x*.  Hence,  x  =  x'  for  some  optimal  x'M 


21 
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fable  2  Job  Moves,  CPU  Times,  and  Final  RELAX  Values  for  Strategies  2,  3,  and  4 


Table  3  Priority  Rule  Heuristic  vs  MCIP  for  Multi-Machine  Reassignment 


